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Abstract 

We consider the quasi-static propagation of a Mode III crack along the interface in a bimaterial plane 
containing a finite array of small line defects (microcracks and rigid line inclusions) . The microdefects are 
arranged to form a channel around the interface that can facilitate (or prevent) the crack propagation. 
pH The two dissimilar elastic materials are assumed to be weakly bonded, so that there is no kinking of 

the main crack from the straight path. On the basis of asymptotic formulae obtained by the authors, 
I the propagation is analysed as a perturbation problem and the incremental crack advance is analytically 

i^h derived at each position of the crack tip along the interface relative to the position of the defects. 

C$ Numerical examples are provided showing potential applications of the proposed approach in the analysis 

of failure of composite materials. Extension to the case of infinite number of defects is discussed. 
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\D 1 Introduction 

m 

OO There is a vast literature devoted to the problem of the interaction between a stationary main crack and 

microdefects, in particular microcracks. The methods of solution range from complex potentials and the 
theory of analytic functions (0 HJ |2]), singular integral equations ([HI HH [15j [17j [16]), and the weight 
function and dipole matrix approach ([TJ |SJ [7])- However, the analysis of crack propagation, even in the 
quasi-static case, requires also the analysis of singular perturbations produced by a small advance of the 

• • crack tip across the microdefect field ( |TUl IT^] ) . 

The symmetric and skew-symmetric weight functions for a two-dimensional interfacial crack, recently 
derived by [10) . can be efficiently employed to analyse the propagation of cracks along the interface in 

$_i heterogeneous materials with small defects. In particular, the notion of skew-symmetric function is essential 

for the problem under consideration , since the "effective" loading produced by the defects on the crack faces 
is not symmetrical, in general. 

We consider the quasi-static propagation of a Mode III crack along the interface in a bimaterial plane 
containing a finite array of small line defects. The microdefects are arranged to form a channel around the 
interface that can facilitate (or prevent) the crack propagation. The defect in question are microcracks and 
rigid line inclusion. The two elastic half-planes are assumed to be weakly bonded, so that there is no kinking 
of the main crack from the straight path. The defects are small compared to the distance from the crack 
tip and the distance between each other and, consequently, the interaction between them can be neglected. 
Using the superposition principle, each defect can be analysed separately and then the solution found by 
summing the contributions for all defects. 
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Fig. [TJshows a main crack at the interface between two dissimilar elastic material interacting with a small 
line defect at a finite distance from the crack tip. The shear moduli for the upper and lower half-planes 
are denoted by /i + and respectively. The line defect, which can be a microcrack or a movable rigid line 
inclusion, has length 2el, where e > is a small dimensionless parameter, and its orientation with respect to 
the £i-axis is determined by the angle a. The centre of the defect is denoted by Y, having polar coordinates 
d and if. The structure is loaded at the crack surfaces by two symmetrical point forces placed at a distance 
a from the crack tip. The crack advance is denoted by e 2 0. 




Figure 1 : A small line defect interacting with an interfacial crack in a bimaterial plane. 



The detailed asymptotic analysis for this problem can be found in |12j . where the approach based on 
weight function and dipole matrix is used. In particular, the dipole fields for different types of defect are 
constructed and used together with symmetric and skew-symmetric weight functions for an interfacial crack 
(|10j) to derive the corresponding perturbation of the stress intensity factor acting at the main crack tip. 

In the next section, we summarise the asymptotic formulae needed for the analysis of propagation of the 
crack in a finite channel of defects, which is presented in Section [3] 



2 Asymptotic formulae 

The propagation of the main crack is described as a perturbation problem, in which the stress intensity 
factor is expanded as follows 
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where = ^2/ira F is the stress intensity factor for the unperturbed crack, AKf u is the singular pertur- 
bation due to a small advance e 2 <p of the crack tip along the interface, and Tl^—i ^^ni ^ s ^ ne perturbation 
produced by a finite number TV of defects. 

Following |12| . the asymptotic formulae for the perturbation terms in (|T|) are given by 
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where Ajjj = —^2/na F/a is the coefficient in the second-order term asymptotics of the unperturbed 
solution, Vu(°\Yj) is the gradient of the unperturbed solution evaluated at Yj, the centre of the defect, 
and Adj is the dipole matrix associated with the type of defect: 
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for a microcrack: 



Mi 



for a movable rigid line inclusion : 
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The full-field solution uq for a Mode III crack in a bimaterial plane without defects can be found in [7] . 

Assuming that the crack is propagating quasi-statically along the interface, the stress intensity factor 
remains constant and equal to the critical value, so that the perturbation term in parenthesis in ([TJ equals 
zero and we get the formula 
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T, AK m> ( 5 ) 
Mil j=i 

which gives the incremental crack advance at each position of the crack tip, driven by the presence of the 
defects, <fi > 0. The procedure is iterated and the total crack elongation at arrest (corresponding to the first 
occurrence of the condition <p < 0) is computed as 



x{K) 
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where K is the number of iterations. 



3 Numerical computations 

We consider the propagation of the main crack driven by a finite number of small line defects arranged to 
form a channel around the interface. Two different examples are discussed below, one consisting of only 
microcracks and the second consisting of microcracks and rigid line inclusions. 

The first example is shown in Fig. [2^, upper part, consisting in 18 microcracks arranged in a 2 x 9 array. 
All the microcracks have the same length, 2el = 0.2, and distance from the interface, h = 1.2, measured from 
the centre of the microcracks. The parameter a denotes the inclination of the microcracks in the first row, 
while microcracks in the second row have inclination a — 7r/2, so that microcracks in the lower half-plane are 
perpendicular to the microcracks in the upper half-plane. The microcracks, above and below the interface, 
are separated one from the another by a constant distance, w = 1. The crack is loaded by a two-point 
symmetrical force applied at a distance a — 0.5 from the crack tip. We analyse first the case of a weak 
interface in a homogeneous plane, that is the contrast parameter r\ = — + is set equal to 

zero. 

In the middle part of the figure we show the corresponding shielding-amplification diagram. We have 
shielding effect when the defects produce a decrease of the stress intensity factor, J2j Ai^niV-R'm < 0, or 
amplification effect in the opposite case, 

J2j AK$/K$ > (E0G3IS]). On the horizontal axis we specify 
the position x of the crack tip along the interface, measured from the point of application of the concentrated 
forces, while the vertical axis stands for the angle a (0 < a < tt). For three values of a, namely a — tt/4, 
7t/2, 37r/4, we report in the lower part of the figure the incremental crack advance tf> as a function of the 
crack tip position x. In the initial position, x — 0.5, the perturbation of the stress intensity factor is positive 
(amplification) and consequently the crack propagates, driven by the defects, and it stops when a neutral 
configuration is reached (line of separation between amplification and shielding regions) . It appears that it 
is not possible to propagate the crack up to the end of the channel. Results show that the value a = 3tt/4 
is the most favourable for the crack propagation. 
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(a) microcracks 



(b) rigid line inclusions/microcracks 



Figure 2: Upper part: (a) Main crack propagating along a weak interface between identical materials containing a 2x9 
array of microcracks. (b) Main crack propagating along a weak interface between identical materials containing a 2x9 array 
consisting of nine rigid line inclusions in the upper half-plane and nine microcracks in the lower half-plane. Middle part: 
Shielding-amplification diagram in the a vs. x plane. Lower part: Incremental crack advance <f> as a function of the crack tip 
position x for three values of a. 



The geometric configuration for the second example is shown in Fig. [2]d, upper part. It consists of a row 
of nine rigid line inclusions in the upper half-plane and a row of nine microcracks in the lower half-plane. 
All defects have the same length 2el = 0.2, inclination a and distance from the interface h = 1.2 and they 
are horizontally separated one from the another by the same distance w — 1. The main crack is loaded by 
a two-point force applied at a distance a = 0.5 from the crack tip. In this case, the shielding-amplification 
diagram shows no dependence on the orientation of the defects, angle a. Consequently, the total crack 
elongation before arrest is independent of a and the crack stops approximately in the middle of the array, 
see the three propagation curves in the lower part of the figure. 

The effect of the position of the two-point force with respect to the crack tip is illustrated in Fig. [3] Fig. 
[3^ corresponds to the case shown in Fig. |2^i but with two-point force applied at distance a = 100 from the 
crack tip, instead of a = 0.5. Fig. [HJd corresponds to the case shown in Fig. but with two-point force 
applied at distance a — 100 from the crack tip. 

Two interesting observations can be made. First, we may note that the "speed" of the crack advance is 
much higher compared to the case when the force is applied close to the crack tip (the incremental crack 
advance (f> is higher by two orders). This can be explained by noticing that, from the unperturbed solution 
uq reported in [7], we deduce 

W°> = 0(a- 1 ' 2 ), <j } = 0(a- 3 / 2 ), a^oo, (7) 

and consequently 

cf) = 0(a), a—too. (8) 
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(a) microcracks (b) rigid line inclusions/microcracks 



Figure 3: Loading applied at large distance from the crack tip. Figure (a) corresponds to the case shown in Fig.|2a but with 
two-point forces applied at distance a = 100 from the crack tip. Figure (b) corresponds to the case shown in Fig. [2p but with 
two-point forces applied at distance a = 100 from the crack tip. 



Second, the crack propagates up to the centre of the channel, x = 5, in the first case, see Fig. [3^,, while it 
propagates across and beyond the whole array in the second case, see Fig.[3]3. Furthermore, the propagation 
curves in Fig. [3]3 show perfect symmetry with respect to the centre of the array. 

All these observations can be confirmed analytically by deriving the asymptotic approximation of <p 
as a — > oo. In particular, we assume the configuration shown in Fig. [4] The crack tip is situated in 
correspondence of the line joining two defects on opposite sides of the interface. Moreover, we number the 
defects with index j running from 1 to N + for the defects ahead of the crack tip, and from — N~ to —1 for 
the defects behind the crack tip. The defects right above and below the crack tip are numbered j = 0. 
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Figure 4: Finite array of defects horizontally separated by a length w from each other and by a length h from the interface. 
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For the case of Fig. [3^, microcracks above and below the interface, we use the formula valid for a single 
microcrack situated at the point Y± (derived in \l'2\ ) 
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By summation over all the defects and assuming that the materials are identical /i + = we obtain as 

a — y oo 
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(10) 



It is noted that the contribution of the defects behind the crack tip, second sum in ( 10 ), equals, but with the 



opposite sign, the contribution of defects ahead of the crack tip, first sum. We deduce that the incremental 
crack advance <j) becomes zero when the crack tip reaches the centre of the array, as confirmed by the 
numerical results in Figj3^. 

For the case of Fig. [3p^ rigid line inclusions above and microcrack below the interface, we use again the 
formula ^ and also the formula valid for a single rigid line inclusion situated at the point Y± (derived in 

AKttt 1 I 2 W=c . /3<p 
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thus obtaining, for identical materials, as a 
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In this case, the contribution of the defects behind the crack tip comes with the same sign as the contribution 
of the defects ahead of the crack tip, which explains the perfect symmetry about the centre of the array 



observed in Fig. |3p. Moreover, since cos 2a is a common factor in (12 1, this formula also explains the constant 
thickness of the amplification region, tt/A < a < 3n/ 4, in Fig. [3J). 

Furthermore, it is possible to extend the analysis to the case of an infinite number of defects ahead of 



the crack tip, N + — s- oo, as the corresponding series in (10) and (12 1 are convergent. In the latter case, we 
obtain an explicit formula as a 
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Finally, we would like to discuss the effects of inhomogeneity on the crack propagation. For a bimaterial 
plane and microcracks arranged as shown in Fig. |2^i, results (not reported here) show that the shielding- 
amplification diagram and the propagation curves are very similar to the case of a homogeneous plane, 
reported in Fig. [2^. This result is expected, since the loading is symmetrical and the defects configuration is 
also symmetrical (apart from the inclination of defects, which is symmetrical only for the angle a = 37r/4). 

We consider now the case of rigid line inclusions above and microcracks below the interface, as shown in 
Fig.|2]D, however now for a bimaterial plane with contrast parameter r\ = —0.67. In this case results, shown 
in Fig. [5j are completely different from the case of a homogeneous plane, compare with Fig. [2]d. Now the 
shielding-amplification diagram shows dependence on the angle a of orientation of the defects with respect 







to the si-axis. Consequently, the total crack elongation at arrest is different for different a's, as shown by 
the three propagation curves in the lower part of the figure. Also, if we swap the two materials, such that 
i] = 0.67, results (not reported here) show that the crack cannot propagate, since a shielding effect takes 
place instead of amplification effect at the initial position of the crack tip. 
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Figure 5: Effects of the inhomogeneity. The figure corresponds to the case shown in Fig. [2J3 but for a bimaterial plane with 
contrast parameter r) = —0.67. 



4 Conclusions 

The problem of propagation of a main interfacial crack in a bimaterial plane containing a finite channel of 
small line defects and subject to antiplane loading conditions is analysed on the basis of asymptotic formulae 
derived by the authors. Under the assumption of quasi-static propagation, the shielding-amplification effect 
of the microdefects on the main crack tip is obtained at each position of the crack tip and the incremental 
crack advance is derived asymptotically. 

Two numerical examples are shown for two different microdefect arrangements. Results show that the 
channel of microdefects can facilitate (amplification) or prevent (shielding) the propagation of the main 
crack. In particular, the orientation of the microdefect can be tuned such that the main crack stops at a 
given position along the channel. We discuss the possibility to extend the asymptotic formulation to the 
case of a channel with infinite number of small defects. 
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